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Running title: A novel MERS-related CoV 


Summary: The discovery of Hp-BatCoV HKU25 bridges the evolutionary gap between MERS-CoV 
and existing bat viruses, and suggests that bat viruses may have evolved to generate MERS-CoV 


through modulation of the spike protein for binding to hDPP4. 


Downloaded from https://academic.oup.com/jid/advance-article-abstract/doi/10.1093/infdis/jiy018/4810771 


by guest 
on 20 January 2018 


Abstract 


Although bats are known to harbor MERS-CoV-related viruses, the role of bats in the 
evolutionary origin and pathway remains obscure. We identified a novel MERS-CoV-related 
betacoronavirus, Hp-BatCoV HKU25, from Chinese pipistrelle bats. While being closely related to 
MERS-CoV in most genome regions, its spike protein occupies a phylogenetic position between 
that of Ty-BatCoV HKUA and Pi-BatCoV HKUS. Since Ty-BatCoV HKU4 but not Pi-BatCoV HKU5 
can utilize MERS-CoV receptor, hDPPA, for cell entry, we tested the ability of Hp-BatCoV HKU25 
to bind and utilize hDPP4. HKU25-RBD can bind to hDPP4 protein and hDPP4-expressing cells, 
but with lower efficiency than that of MERS-RBD. Pseudovirus assays showed that HKU25-spike 
can utilize hDPP4 for entry to hDPP4-expressing cells, though with lower efficiency than that of 
MERS-spike and HKU4-spike. Our findings support a bat origin of MERS-CoV and suggest that bat 


coronavirus spike proteins may have evolved in a stepwise manner for binding to hDPP4. 


Keywords: Middle East Respiratory Syndrome Coronavirus, Spike glycoprotein, Dipeptidyl 


peptidase 4, Hypsugo bat 
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Introduction 


The Middle East Respiratorv Svndrome (MERS) has affected 27 countries in four 
continents with 2090 cases and a fatalitv rate of 34.996 since its emergence in 2012. The 
etiological agent, MERS coronavirus (MERS-CoV), belongs to Betacoronavirus lineage C (1, 2] and 
utilizes human dipeptidyl peptidase 4 (hDPP4) as receptor for cell entry [3]. While dromedaries 
are likely the immediate animal source of the epidemic [4-6], bats also harbor MERS-CoV-related 
viruses which may suggest a possible bat origin [7-13]. However, the evolutionary pathway and 
direct ancestor of MERS-CoV remains obscure. In particular, there is an evolutionary gap 


between MERS-CoV and related bat viruses. 


Since the SARS epidemic, numerous novel CoVs have been discovered [14-16], with bats 
uncovered as an important reservoir for alphacoronaviruses and betacoronaviruses [17-21]. 
When MERS-CoV was first discovered, it was most closely related to Tylonycteris bat CoV HKU4 
(Ty-BatCoV HKUA) and Pipistrellus bat CoV HKU5 (Pi-BatCoV HKU5) previously discovered in 
Lesser bamboo bat (Tylonycteris pachypus) and Japanese pipistrelle (Pipistrellus abramus) 
respectively in Hong Kong [1, 7-10, 22]. The spike of Ty-BatCoV HKUA, but not that of Pi-BatCoV 
HKU5, was able to utilize the MERS-CoV receptor, hDPP4 or CD26, for cell entry (3, 23]. 
Subsequently, three other lineage C betacoronaviruses, Coronavirus Neoromicia/PML- 
PHE1/RSA/2011 (NeoCoV), BtVs-BetaCoV/SC2013 and BatCoV PREDICT/PDF-2180 were also 
detected in vesper bats from China or Africa [11-13, 24]. A lineage C betacoronavirus, Erinaceus 
CoV VMC/DEU, has also been found in European hedgehogs [25]. This is interesting because 
hedgehogs are phylogenetically closely related to bats. MERS-CoV can infect bat cell lines and 
Jamaican fruit bats [25, 26], further suggesting that bats may be the primary host of MERS-CoV 


ancestors. 
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Although NeoCoV represents the closest bat counterpart of MERS-CoV in most genome 
regions, its spike (S) protein is genetically divergent from that of MERS-CoV [11], suggesting an 
evolutionary gap between existing MERS-CoV and bat viruses and an immediate ancestor of 
MERS-CoV yet to be discovered. To identify the potential bat origin and understand the 
evolutionary path of MERS-CoV, we collected bat samples from various regions in China. Diverse 
CoVs were detected, including a novel lineage C betacoronavirus from Chinese pipistrelle 
(Hypsugo pulveratus), which can utilize hDPP4 for cell entry. The results support a bat origin of 


MERS-CoV and suggested stepwise evolution of spike protein in hDPP4 binding. 
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Materials and methods 


Ethics statement. Bat samples were collected bv Guangdong Institute of Applied 
Biological Resources, Guangzhou, China, in accordance with guidelines of Regulations for 
Administration of Laboratorv Animals under a license from Guangdong Entomological Institute 


Administrative Panel on Laboratorv Animal Care. 


Detection of CoVs from bats. Samples were collected from bats captured from various 
locations in seven provinces of China (Figure 1) during 2013-2015 using procedures described 
previously (27, 28]. Viral RNA extraction was performed using QIAamp Viral RNA Mini Kit 
(QIAgen, Hilden, Germany). CoV detection was performed by Reverse-transcription polymerase 
chain reaction (RT-PCR) targeting a 440-bp fragment of RNA-dependent RNA polymerase (RdRp) 
gene using conserved primers (5'-GGTTGGGACTATCCTAAGTGTGA-3' and 5'- 
ACCATCATCNGANARDATCATNA-3') as described previously [16]. A phylogenetic tree was 


constructed with maximum likelihood method using GTR+G+I substitution model by MEGA 6.0. 


Viral culture. The two Hp-BatCoV HKU25 samples were subject to virus isolation in Vero 
E6 (ATCC CRL-1586), Huh-7 (JCRBO403), PK15 (ATCC CCL-33) and Rousettus lechenaultii primary 


kidney cells (in-house) as described previously [29]. 


Complete genome sequencing and analysis of Hp-BatCoV HKU25. Two Hp-BatCoV 
HKU25 complete genomes were sequenced according to our published strategy [27]. A total of 
75 sets of primers, available on request, were used for PCR. The assembled genome sequences 
were compared to those of other CoVs using the comprehensive coronavirus database CoVDB 


(http://covdb.microbiologv.hku.hk) [30]. The time of the most recent common ancestor 
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(tMRCA) was estimated based on ORF1ab sequences, using uncorrelated exponential distributed 


relaxed clock (UCED) model in BEAST version 1.8 (http://evolve.zoo.ox.ac.uk/beast/) [31]. 


Cloning of recombinant S1-receptor-binding-domain (RBD) proteins. The S1-RBD 
sequences of of Hp-BatCoV HKU25 (residues 374-604) and MERS-CoV (residues 367-606) were 
cloned into mammalian expression vector pCAGGS containing signal peptide (CD5) and C- 
terminal Fc tag from mouse IgG2a (mFc) (32, 33]. The expression plasmids were transiently 
transfected into human embryonic kidney HEK293T cells (ATCC CRL-3216). The recombinant 
HKU25-RBD-mFc and MERS-RBD-mFc proteins were purified by protein A-based affinity 


chromatography. 


Protein binding with flow cytometry and fluorescence-activated cell sorter (FACS) 
analysis. Huh7 (normal or DPP4 knockdown using small interfering RNA (siRNA)) or 293T 
(normal or transfected with DPP4-expressing plasmid) cells were incubated with 10 ug/ml 
MERS-RBD-mFc or 40 pg/ml HKU25-RBD-mFc at 4°C for 1 h. Cells were then stained with Alexa 
Fluor 488-conjugated goat anti-mouse IgG on ice for 30 min. Protein-to-cell binding was 


analyzed using BD FACS LSRII instrument (BD Bioscience, East Rutherford, New Jersey, USA). 


Immunostaining and confocal microscopy. Huh7 cells were fixed on glass coverslips and 
incubated with 50 ug/ml HKU25-RBD-mFc or 20 ug/ml MERS-RBD-mFc in PBS at 4'C for 1 h, 
followed by staining with Alexa Fluor 488-conjugated goat anti-mouse or anti-rabbit IgG. Cell 
nuclei were stained using 4',6-diamidino-2-phenylindole (DAPI) in mounting medium. Images 
were acquired with 63x oil objectives using a Zeiss LSM510 Meta laser scanning confocal 


microscope. 


Knockdown of hDPP4 expression using siRNAs. siRNA duplexes against hDPP4 (5'- 


UGACAUGCCUCAGUUGUAUU-3’) were synthesized by Nucleic Acids Center at National Institute 
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of Biological Sciences, Beijing, China, with non-targeting siRNA as negative control (Ctrl-si). Ten 
picomoles of siRNA were transfected into Huh7 cells with Lipofectamine RNAiMax (Invitrogen). 
Knockdown efficiency was determined by quantitative qRT-PCR analysis using primers specific 
for hDPP4 (5'-CCTGCTTCTATGTTGATA -3'; 5'-CGAATAGTTCTGAATCCT -3') and western blot 
analysis using anti-hDPP4 antibody (Abcam, Cambridge, United Kingdom). The mRNA levels 

of target genes were normalized to that of glyceraldehyde 3-phosphate dehydrogenase (gapdh) 


gene [34]. 


Immunoprecipitation. To identify the direct interaction between MERS-RBD-mFc or 
HKU25-RBD-mFc and hDPP4, HEK 293T cells were transfected with hDPP4-expressing plasmids 
and lysed with RIPA buffer containing 1x protease inhibitor cocktail (Roche) 48 h after 
transfection. Cell lysates were incubated with purified MERS-RBD-mFc or HKU25-RBD-mFc and 
Dynal protein A Sepharose beads at 4*C overnight. The bound fractions of immunoprecipitates 
(IP) and total cell lysate (as input) were analyzed by western blot with anti-mFc, anti-hDPP4 or 


anti-GAPDH antibodies. 


Pseudovirus production. Retroviruses pseudotyped with MERS-CoV, Ty-BatCoV HKUA, 
Pi-BatCoV HKU5 and Hy-BatCoV HKU25 S proteins were packaged by HEK293FT cells (R70007, 
Invitrogen). Briefly, plasmid containing the respective CoV S gene was co-transfected with a 
plasmid containing luciferase gene but env-defective HIV-1 (pNL 4-3.Luc.RE) into 293FT cells 
using Lipofectamine 2000 (Invitrogen). Culture supernatant was concentrated with 5x PEG-it 
virus precipitation solution (SBI). For mock pseudoviruses (Aenv) bearing no S protein, empty 


plasmid was co-transfected with pNL 4-3.Luc.RE. 


Pseudovirus cell entry assay. HEK293T cells were transfected with plasmid containing 


hDPP4 gene and empty plasmid (as mock-transfected control) by Lipofectamine 2000. 
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Pseudoviruses bearing CoV S proteins were treated by 100ug/ml Tosyl phenylalanyl 
chloromethyl ketone (TPCK)-treated trypsin at 37°C for 30 min prior to infection. After trypsin 
inactivation, pseudovirus infections were performed by spinning at 1200 g at 4°C for 2 h and 
incubation at 37*C for 5 h. Cells were then incubated for 72 h and Ivsed for luciferase activity 
determination using Luciferase Assay System (Promega, Fitchburg, USA). To test for inhibition of 
pseudovirus-mediated cell entry by anti-hDPP4 antibodies, HEK293T cells transfected with 
hDPP4 were pre-incubated with 10 ug/ml anti-hDPP4 polyclonal antibodies (R&D systems) at 


37*C for 1 h before pseudovirus infection. 


Structural modelling of Hp-BatCoV HKU25 RBD. The model of HKU25-RBD and HKU5- 
RBD was built with the crystal structure of MERS-RBD/hDPP4 using SWISS-MODEL with default 
parameters and analyzed using Discovery Studio visualizer (Accelrys, San Diego, USA), and the 
Ramachandran plot were examined to ensure that the structure of the models were not in any 
unfavorable region. The models of HKU4-RBD and HKU5-RBD were also built as positive and 
negative controls respectively with the same parameters, and were superimposed for 


comparison. 


Nucleotide sequence accession numbers. The nt and genome sequences of CoVs 
detected in this study have been lodged within GenBank under accession no. KX442564, 


KX442565, and KX447541 to KX447565. 
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Results 


Detection of CoVs in bats and discoverv of a novel lineage C betacoronavirus from Chinese 


pipistrelle. 


A total of 1964 alimentarv samples from bats belonging to 19 different genera and 44 
species were obtained from seven provinces of China. RT-PCR for a 440-bp fragment of RdRp 
gene of CoVs was positive in samples from 29 bats of five species belonging to four genera 
(Figure 1 and Supplementary Table 1). Sequence analysis showed that four samples contained 
alphacoronaviruses, five contained lineage B betacoronaviruses and 20 contained lineage C 


betacoronaviruses (Supplementary Figure 1). 


Of the 20 lineage C betacoronavirus sequences, 18 sequences from Tylonycteris 
pachypus possessed 9696 nt identities to Ty-BatCoV HKUA. The other two lineage C 
betacoronavirus sequences (YD131305 and NL140462) showed <86% nt identities to MERS-CoV 
or other lineage C betacoronaviruses, suggesting a potentially novel lineage C betacoronavirus 
closely related to MERS-CoV (Supplementary Table 1 and Supplementary Figure 1). Both 
samples were collected from Chinese pipistrelle (Hypsugo pulveratus) bats, belongs to the family 
Vespertilionidae, captured in Guangdong Province (Figure 1). We proposed this novel CoV to be 
named Hypsugo pulveratus bat coronavirus HKU25 (Hp-BatCoV HKU25). Attempts to passage 


Hp-BatCoV HKU25 YD131305 and NL140462 in cell cultures were not successful. 


Genome features of Hp-BatCoV HKU25. 


The complete genome sequences of the two Hp-BatCoV HKU25 strains, YD131305 and 
NL140462, were determined, with genome features similar to MERS-CoV including conserved 


ORF4a and ORF4b (Supplementary Table 2, Supplementary Table 4, Supplementary Figure 2 and 
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Supplementarv Figure 3). Thev shared 95.996 overall nt identities, while possessing 82.096, 73.2- 
73.996, 73.596 and 69.396 nt identities to the genomes of BtVs-BetaCoV/SC2013, human/camel 
MERS-CoVs, NeoCoV and Ty-BatCoV HKUA respectively. Comparison of the seven conserved 
replicase domains for CoV species demarcation showed that Hp-BatCoV HKU25 represents a 
novel species under Betacoronavirus lineage C (Supplementary Table 5), with the concatenated 
sequence being most closely related to that of BtVs-BetaCoV/SC2013 with 88.596 amino acid (aa) 


identities. 


Phylogenetic and molecular clock analysis. 


Phylogenetic trees constructed using RdRp, ORF1, S1 and N sequences of Hp-BatCoV 
HKU25 are shown in Figure 2. Hp-BatCoV HKU25 was most closely related to BtVs- 
BetaCoV/SC2013, forming a distinct branch among lineage C betacoronaviruses. In RdRp, ORF1 
and N genes, MERS-CoVs were most closely related to NeoCoV followed by the branch formed 
by Hp-BatCoV HKU25 and BtVs-BetaCoV/SC2013. In contrast, in S1 region, MERS-CoVs were 
most closely related to Ty-BatCoV HKUA, followed by the branch formed by Hp-BatCoV HKU25 
and BtVs-BetaCoV/SC2013, but was only distantly related to NeoCoV. Hp-BatCoV HKU25 and 
BtVs-BetaCoV/SC2013 thus represent close relatives of MERS-CoV, while they occupied a 


position in between Ty-BatCoV HKUA and Pi-BatCoV HKUA in relation to MERS-CoV in S1 region. 


Using the uncorrelated relaxed clock model on ORF1ab, tMRCA of human and camel 
MERS-CoVs was dated to 2009.56 [Highest Posterior Density (HPD), 2006.8-2011.3], while that 
of MERS-CoV, NeoCoV, Hp-BatCoV HKU25 and BtVs-BetaCoV/SC2013 was dated to 1939.32 


(HPDs, 1899.5-1969.0) (Supplementary Figure 4). 
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Sequence analysis of Hp-BatCoV spike protein. 


MERS-CoV utilizes hDPP4, a type ll transmembrane protein, as receptor for initiation of 
infection [3]. The S1 domain responsible for hDPP4 receptor binding is located in a C-terminal 
240-residue RBD that contains the receptor binding motif (RBM) which engages the receptor 
[35]. Using binding and pseudovirus assays, it was shown that Ty-BatCoV HKUA S, but not Pi- 
BatCoV HKUS S, can bind to and utilize hDPPA for cell entry (23, 36]. Since phylogenetic analysis 
placed Hp-BatCoV HKU25-S1 at a position between Ty-BatCoV HKUA and Pi-BatCoV HKUS in 
relation to MERS-CoV, it would be interesting to know if Hp-BatCoV HKU25 may bind to and 
utilize hDPPA for cell entry. As in other CoVs, Hp-BatCoV HKU25-S is predicted to be a type | 
membrane glycoprotein with two heptad repeats. The predicted Hp-BatCoV HKU25-S1-RBD 
shared 53.596 aa identities to that of MERS-CoV, with two short deletions compared to MERS- 


CoV and Ty-BatCoV HKU4 (Figure 3). 


Previous structural studies have identified 12 critical residues (Y499, L506, D510, E513, 
W535, E536, D537, D539, Y540, R542, W553 and V555) for hDPP4 binding in MERS-CoV (23, 37]. 
In Ty-BatCoV HKUA, five (Y503, L510, E518, E541 and D542 corresponding to V499, L506, E513, 
E536 and D537 in MERS-RBD) residues were conserved, which may allow binding to hDPP4. In 
Pi-BatCoV HKU5, one of the 12 conserved residues (D543 corresponding to D537 in MERS-RBD) 
was found. In Hp-BatCoV HKU25, one residue (R546 in strain YD131305/R547 in strain NL140462 
corresponding to R542 in MERS-RBD) was conserved in both strains and an additional residue 


(V554 corresponding to V555 in MERS-RBD) was conserved in strain NL140462 (Figure 3). 


Hp-BatCoV HKU25-S1-RBD binds hDPP4 but with lower efficiency than MERS-CoV S1-RBD. 


To examine the ability of Hp-BatCoV HKU25-S to bind hDPP4, we expressed and purified 
hDPP4 and the S1-RBD domains of Hp-BatCoV HKU25 (residues 374-604) and MERS-CoV 
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(residues 367-606) using procedures described previously (32, 33]. We first tested the binding 
efficiency of the S1-RBD domains to Huh7 cells (human hepatocellular carcinoma cells with 
endogenous hDPP4 expression) using flow cytometry. HKU25-RBD can bind to Huh7 cells, 
although the observed fluorescence shift was smaller than MERS-RBD (Figure 4, panel A). This 
indicates that HKU25-RBD can bind to hDPP4-expressing Huh7 cells with lower binding efficiency 
than that of MERS-RBD. To confirm that the binding is mediated by hDPP4, we obtained Huh7 
cells with small interfering RNA (siRNA) knockdown of hDPP4 (confirmed by mRNA expression 
and western blot) (Figure 4, panel C). A significant reduction of fluorescence shift was observed 
in both HKU25-RBD- and MERS-RBD-mediated binding to hDPP4-knockdown Huh7 cells when 
compared to hDPP4-expressing Huh7 cells (Figure 4, panel A). Moreover, HKU25-RBD and MERS- 
RBD could only bind to HEK293T cells (lacking endogenous hDPP4 expression) after transfection 
with hDPP4-expressing plasmid, although the binding efficiency to hDPP4-expressing HEK293T 
cells was lower for HKU25-RBD than MERS-RBD (Figure 4, panel B). Second, we also confirmed 
the binding of HKU25-RBD to Huh7 cell surface by confocal microscopy (Figure 4, panel D). Third, 
immunoprecipitation assays showed that hDPP4 protein can be specifically pulled down by both 
MERS-RBD and HKU25-RBD (Figure 5). These results indicated that HKU25-RBD can bind to 


hDPP4 on cell surface, but with lower efficiency than MERS-CoV-RBD. 


HKU25 pseudovirus can utilize hDPP4 for cell entry but with lower infection efficiency than 


MERS and HKU4 pseudoviruses. 


To determine if Hp-BatCoV HKU25-S can mediate viral entry into hDPP4-expressing 
human cells, we performed HKU25-S-mediated pseudovirus entry assay. Since the S protein of 
Ty-BatCoV HKU4 but not that of Pi-BatCoV HKU5 can utilize hDPP4 for cell entry [36], we 


included HKU4-S, HKU5-S and MERS-S mediated pseudovirus entry assays for comparison. 
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Pseudovirus assays were used because isolation of live Hp-BatCoV HKU25 was not successful, as 
with most bat CoVs which are often difficult to culture. Retroviruses pseudotvped with 
luciferase and the respective S proteins were tested for entrv to HEK293T cells with or without 
hDPP4 expression. MERS-S most robustly mediated pseudovirus entry into hDDP4-expressing 
HEK293T cells, followed by HKU4-S and HKU25-S, as shown by luciferase activities measured. All 
three pseudoviruses showed marked increase in luciferase activities in hDDP4-expressing 
HEK293T cells compared to cells without hDPP4 expression (Figure 6). Moreover, anti-hDPP4 
polyclonal antibodies could competitively block HKU25-S, HKU4-S and MERS-S pseudovirus entry 
to hDPP4-expresssing HEK293T cells, further confirming the binding specificity. In contrast, 
HKU5-S and control retroviruses not pseudotyped with S did not mediate pseudovirus entry into 
hDPP4-expressing HEK293T cells (Figure 6). These results showed that hDPP4 is a possible 
functional receptor for Hp-BatCoV HKU25, although cell entry may be less efficient than Ty- 


BatCoV HKU4 and MERS-CoV. 


Structural modelling of RBD-hDPP4 binding interphase 


To predict the RBD-hDPP4 binding-interface, the structures of HKU25-, MERS-, HKU4- 
and HKU5-RBDs were modelled with that of hDPP4 using homology modelling. The seguence 
identity between HKU25-RBD and MERS-RBD (template) was 25096 and the RBD-hDPP4 interface 
for all RBDs was similar (Supplementary Figure 5), except that only MERS-RBD and HKU4-RBD 
possess the extended loop between B6 and B7 involved in interaction with hDPP4 [23]. A 
negative-charge residue, E536, located in MERS-RBD, corresponding to E541 in HKU4-RBD, can 
interact with the carbohydrate moiety of hDPP4, whereas HKU5-RBD contains a positive-charge 


residue, R542, and HKU25-RBD contains an uncharged residue, T540/A541 at the corresponding 
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position. These findings supported that the binding of HKU25-RBD to hDPP4 may be weaker 


than that of MERS-RBD and HKU4-RBD but stronger than that of HKU5-RBD. 
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Discussion 


The novel lineage C betacoronavirus, Hp-BatCoV HKU25, helps to fill the evolutionary 
gap between existing bat viruses and MERS-CoV, and offers new insights into the evolutionary 
origin of MERS-CoV. Hp-BatCoV HKU25 shared similar genome features with MERS-CoV, 
including the conserved ORF4a and ORF4b with predicted domains for dsRNA binding and 
antagonizing interferon signals respectively [38, 39]. Phylogenetically, Hp-BatCoV HKU25, 
together with BtVs-BetaCoV/SC2013, was closely related to MERS-CoV and NeoCoV in most 
genome regions, suggesting that these viruses share a common ancestral origin. While the S1 of 
NeoCoV is only distantly related to MERS-CoV, the S1 of Hp-BatCoV HKU25 was at a 
phylogenetic position closely related to MERS-CoV, only second to Ty-BatCoV HKU4. On the 
other hand, the S1 of NeoCoV is most closely related to Erinaceus CoV from European 
hedgehogs. Since NeoCoV was detected in an African bat, it is more likely a recombinant virus 
between bat and hedgehog CoVs in Africa. Moreover, it was shown that the S of BatCoV 
PREDICT/PDF-2180, which is closely related to NeoCoV in all genome regions, cannot mediate 
entry to hDPP4-expressing cells [13]. This further supported that NeoCoV and PREDICT/PDF- 
2180 are unlikely the immediate ancestors of MERS-CoV. On the other hand, Hp-BatCoV HKU25 
and related viruses may represent close relatives to the immediate ancestor of MERS-CoV, 


based on its phylogenetic position in all genome regions including S protein. 


The ability of Hp-BatCoV HKU25 to utilize hDPP4 for cell entry suggests that the S 
protein of related bat viruses may have evolved to cross the species barrier during the 
emergence of MERS-CoV. Using binding and pseudovirus assays, we demonstrated the ability of 
HKU25-S to bind to and utilize hDPP4 for cell entry, though with infection efficiency lower than 


that of MERS-S and HKU4-S. This is not only in line with the phylogenetic position of HKU25-S1 
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between HKU4-S1 (which can bind to hDPP4) and HKU5-S1 (which cannot), but is also consistent 
with findings from structural modelling. Our results suggested MERS-CoV may have originated 
from bat viruses having acguired a stepwise increasing ability to bind hDPP4 as they evolved. 
Previous molecular dating studies estimated that the time of divergence of MERS-CoVs was 
approximately 2010/2011 [40-43]. The present dating results are in line with such estimation, 
with the tMRCA of MERS-CoVs dated to approximately 2009, and that of MERS-CoV, NeoCoV, 
Hp-BatCoV HKU25 and BtVs-BetaCoV/SC2013 dated to approximately 1939. Therefore, the 
immediate ancestor of MERS-CoV could well have emerged from bats in the last century 


through evolution in its S protein before jumping to camels and humans. 


The evolutionary path of MERS-CoV may be different from that of SARS-CoV. For SARS- 
CoV, the overlapping habitat and geographical distribution of different horseshoe bats in China 
is believed to have fostered viral recombination leading to the epidemic. SARS-CoV is most likely 
a recombinant virus arising from ancestral viruses in horseshoe bats before it jumped to civet 
and then humans [44-47]. In contrast, there is currently no evidence to suggest that MERS-CoV 
is a recombinant virus. A previous study suggested that the genetically divergent S1 in NeoCoV 
may indicate intraspike recombination events involved in the emergence of MERS-CoV [11]. As 
explained above, NeoCoV, rather than MERS-CoV, is more likely a recombinant virus. On the 
other hand, a stepwise evolution of the S protein in gaining the ability to utilize camel and 
human DPP4 may be an important mechanism for interspecies transmission during the 


emergence of MERS-CoV. 


Our results further support a possible bat origin of MERS-CoV and suggest that 
continuous surveillance of bats in the Middle East, Africa and other regions may reveal the 


immediate animal origin of MERS-CoV. The application of similar state-of-the-art molecular 
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studies on naturallv evolving ancestral and intermediate viruses along the evolutionarv path 
mav provide further clues in understanding the mechanism of interspecies transmission of 
emerging viruses, while obviating the risks of generating dangerous mutants using the 


controversial, gain-of-function studies. 
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LEGENDS TO FIGURES 


Figure 1 Map showing various sampling locations in seven provinces of China (Guangxi, Guangdong, 
Shanxi, Zhejiang, Yunnan, Hainan and Guizhou). Sampling locations with Hp-BatCoV HKU25 and other 


CoVs detected are in blue and red respectively. 


Figure 2 Phylogenetic analyses of RdRp, ORF1, S1 and N nucleotide sequences of Hp-BatCoV HKU25 and 
other lineage C betacoronaviruses (B). The trees were constructed by maximum likelihood method using 
GTR-G substitution models respectively and bootstrap values calculated from 1000 trees. Trees were 
rooted using corresponding sequences of HCoV HKU1 (GenBank accession number NC 006577). Only 
bootstrap values >70% are shown. (A) 2775 nt (B) 20694 nt (C) 3740 nt (D) 1167 nt positions respectively 
were included in the analyses. The scale bars represent (A) 20 (B) 20 (C) 10 (D) 10 substitutions per site 
respectively. Human and camel MERS-CoVs are in purple, Ty-BatCoV HKUA and Pi-BatCoV HKUS are in 
blue, NeoCoV is in green and BetaCoV/SC2013 is in pink. The two Hp-BatCoV HKU25 strains, YD131305 
and NL140462, detected in this study are in red and bolded. Structural modeling of the receptor-binding 
domain (RBD) of the spike protein of Hp-BatCoV HKU25 (B). The models of RBDs of HKU4-4S (green), 
HKU25-YD131305 (magenta), HKU25-NL140462 (gold) and HKU5-27S (gray) are shown with hDPP4 
structure (violet) in ribbon diagram. The interface of different RBDs and hDPP4 are zoomed and the 
residues that highlighted in multiple seguence alignment from Fig. 2 are shown in ball-and-stick format 
colored by element (carbon, gray; nitrogen, blue; oxygen, red). Strands of B6 and B7 are labeled in the 


structure of HKU4-4S only. The figures were produced using Discovery Studio visualizer (Accelrys). 


Figure 3 Multiple alignment of the amino acid seguences of the receptor-binding domain (RBD) of the 
spike protein of MERS-CoV and corresponding seguences in Hp-BatCoV HKU25 and other lineage C 


betacoronaviruses. Asterisks indicate positions with fully conserved residues. The two amino acid 
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deletions in Hp-BatCoV HKU25 compared to MERS-CoV and Ty-BatCoV HKUA are indicated with red 
boxes. The 12 critical residues for receptor binding in MERS-CoV are highlighted in different colors. The 
10 residues marked below the alignment are based on (Wang, 2013)[37] and the other two residues 
marked above the alignment are based on (Wang, 2014)[23]. Y499, highlighted in blue, formed 
hydrogen bond with DPP4 residue. L506, W553 and V555, highlighted in green, formed a hydrophobic 
core surrounded by hydrophilic residues D510, E513 and Y540, which are highlighted in yellow. D510 
and E513 also contributed to salt bridge interaction and hydrogen bonding with DPP4 residues. E536, 
D537 and D539, highlighted in pink, formed negative-charged surface. W535, highlighted in grey and 


R542, highlighted in purple, are residues that have strong polar contact with DPP4 residues. 


Figure 4 Binding of HKU25-RBD with human cells was mediated by interacting with hDPP4 receptor. (A) 
FACS analysis of MERS-RBD-mFc (10 ug/ml) and HKU25-RBD-mFc (40 ug/ml) binding to Huh7 cells and 
hDPP4-knockdown Huh7 cells. (B) FACS analysis of MERS-RBD-mFc and HKU25-RBD-mFc binding to 293T 
cells and 293T cells transfected. with hDPP4-expressing plasmid. The shaded area represents the 
secondary antibody control. (C) Determination of siRNA efficiency by qRT-PCR and western blot analysis 
using primers and antibody specific for hDPP4. (D) MERS-RBD-mFc and HKU25-RBD-mFc binding to a 
molecule(s) located on the Huh7 cell surface. MERS-RBD-mFc and HKU25-RBD-mFc were detected by an 
Alexa Fluor 488-conjugated goat anti-mFc antibody. Empty expressing plasmid was used as a negative 


control. 


Figure 5 MERS-RBD-mFc (A) and HKU25-RBD-mFc (B) proteins directly bind with hDPP4. HEK 293T cells 
were transfected with hDPP4-expressing plasmids, MERS-RBD-mFc (A) and HKU25-RBD-mFc (B) proteins 


were used for immunoprecipitation of lysates of HEK 293T cells transfected with hDPP4-expressing or 
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empty plasmids. Empty plasmid was mock-transfected as negative control. hDPP4 was coprecipitated 


from the lysates, as detected by antibody specific for hDPP4. GAPDH was used as a loading control. 


Figure 6 HEK293T cells transfected with empty plasmid or hDPP4 were infected by retroviruses 
pseudotyped with MERS-CoV, Ty-BatCoV HKUA, Pi-BatCoV HKU5 and Hp-BatCoV HKU25 S proteins with 
mock pseudovirus (Aenv) as control. The cells were also preincubated with anti-hDPP4 antibodies to 
test for cell entry inhibition. Cell entry efficiencies were assayed by luciferase activity measurement 


after 72 h. 
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Figure 1. 
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Figure 2. 
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Figure 3. 
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Figure 4. 
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Figure 6. 
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